Transition Voltage Spectroscopy and the nature of vacuum tunneling 
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Transition Voltage Spectroscopy (TVS) has been proposed as a tool to analyze charge transport 
through molecular junctions. We extend TVS to Au- vacuum- Au junctions and study the distance 
dependence of the transition voltage K(d) for clean electrodes in cryogenic vacuum. On the one 
hand, this allows us to provide an important reference for Vt (d)-measurements on molecular junc- 
tions. On the other hand, we show that TVS forms a simple and powerful test for vacuum tunneling 
models. 
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Molecular electronics offers the exciting prospect of 
tuning the conductance of nanometer-sized junctions 
with the tools of synthetic chemistry. In principle, it is 
possible to fabricate small switchable devices from func- 
tional molecules However, the electronic properties 
of a molecule change dramatically when it is contacted by 
metal electrodes, as the coupling significantly shifts and 
broadens the molecular levels Hence, it is of utmost 
importance to understand the influence of the electrodes 
on the molecular levels. The level structure can in prin- 
ciple be studied by measuring the current-voltage (I-V) 
characteristics of molecular junctions; as soon as a molec- 
ular level becomes accessible by the bias window, a step 
in the current can be observed. However, such exper- 
iments are complicated by the fact that the position of 
the molecular level is often a few eV away from the Fermi 
energy Ep. As a consequence, very large electric fields 
need to be applied to probe these levels, often resulting 
in a breakdown of the junction before the first molecular 
level can be accessed. 

Recently, Beebe et al. have presented a new method 
to characterize the electronic structure of molecular junc- 
tions which requires lower electric fields 0, [H| • By plot- 
ting their I-V data in a Fowler-Nordheim (FN) repre- 
sentation, i.e. plotting ln(//y^) versus 1/V [g, they 
observed a minimum at a characteristic voltage Vt. In- 
terestingly, for a broad set of molecules 14 was found to 
be proportional to the position of the highest occupied 
molecular orbital [eh) with respect to Ep, as confirmed 
by ultraviolet photoelectron spectroscopy (UPS). Hence, 
TVS was proposed as a new, simple and powerful spec- 
troscopic tool to study the positions of molecular levels. 

However, the interpretation of TVS is still under de- 
bate f?'-!^] . In the original picture of Beebe et al. a molec- 
ular junction was treated as a simple tunnel junction. 
They interpreted Vt as the potential for which the tunnel 



barrier, tilted by the bias voltage, changes from trape- 
zoidal to triangular. In this picture, Vt would be the same 
for all electrode separations d. In contrast, Huisman et 
al. considered the tunnel barrier model in more detail 
and calculated that Vt should be strongly dependent on 
d 0, Ell • They obtained the following approximation 
for Vt in the case of a trapezoidal vacuum barrier with 
height (m is the electron mass): 
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Remarkably, such a variation of Vt with d is much 
stronger than the weak distance dependence reported for 
junctions containing aliphatic and conjugated molecules 
of different lengths 5] . This statement remains valid even 
when including image forces which modify the potential 
barrier Q. Therefore, Huisman et al. suggested that the 
dependence Vt(d) might form a new tool to distinguish 
molecular junctions from vacuum junctions. However, 
TVS studies on vacuum junctions have not been reported 
thus far. 

To test this proposition we present extensive TVS mea- 
surements on metal-vacuum-metal junctions. The possi- 
bility to accurately vary the tunnel gap between the elec- 
trodes allows us to reveal the basic properties of TVS. We 
show that for Au junctions Vt only has a weak variation 
with electrode spacing d, comparable to the dependence 
reported for molecular junctions. Hence, our results in- 
dicate that the Vt(d) dependence itself cannot be used to 
distinguish molecular junctions from vacuum junctions. 
However, the absolute values of Vt are significantly higher 
for vacuum barriers than for junctions with conjugated 
molecules. Interestingly, our study demonstrates that 
TVS can also be used as a characterization tool for vac- 
uum barriers. The dependence <f)(d) in standard models 
for vacuum tunneling, such as those leading to Eq. [TJ 
can be tested by TVS measurements. 

For a full characterization of TVS on a metal- vacuum- 
metal junction, voltages up to 3 V are to be applied over 
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vacuum gaps as small as «0.3 nm. Such high electric 
fields (and high field gradients) may cause instabilities 
in the tunnel junctions. Hence, junctions are needed 
which are stable in time and kept in a clean environ- 
ment. For this reason we used notched- wire mechani- 
cally controllable break junctions in cryogenic vacuum 
(Tft! 5 K) The electrodes are made of gold, the 

archetypical electrode metal for molecular junctions. In 
addition, the junctions were first optimized by a " train- 
ing" procedure,__Le. by repeatedly opening and closing 
the electrodes [ifi]. We expect that this organizes the 
apex atoms into their strongest bond configuration and 
enhances their stability in high electric fields. The high 
stability and repeatability of the conductance evolution 
is illustrated by the tunnel curve in Fig. [1^. Upon clos- 
ing, the tunnel current increases exponentially until the 
electrodes snap to contact [l^ . Note that the con- 
ductance jumps to a value close to 1 Go (IGo — 26^ /h), 
indicating a clean single-atom contact. After the training 
procedure, the electrodes are separated such that a vac- 
uum gap is created with a zero bias conductance of ssO.Ol 
Gq. This is the starting point for the TVS measurements. 
Subsequently, the vacuum gap is increased stepwise, and 
an I-V curve over the range ±2-3 V is recorded for each 
position. 

A typical example of such an I-V curve is plotted in 
Fig. [TJd. Clearly, the current displays a transition from 
a linear dependence at low voltages to a strongly non- 
linear behavior for voltages > l.bV . As stated earlier, 
this transition can be quantified by scaling the data in a 
Fowler-Nordheim representation. This is shown in Fig. 
[T];, here for positive bias voltage only. Let us first discuss 
the upper curve. This curve is measured for a small tun- 
nel gap with a zero bias conductance of «0.02 Go. It has 
a well-defined minimum at 1.64 V that determines Vt. In 
total, 34 curves are plotted, corresponding to 34 differ- 
ent electrode separations (equally spaced by RiO.02 nm). 
When we increase the electrode separation, a shift in Vt 
can be observed. The transition voltage first increases 
with distance and is at a maximum after stretching by 
nm (fifth curve). Vt has now increased to 1.9 V. For 
even larger gaps Vt decreases again to a value of 1.55 V 
after stretching by wO.G nm. 

Before we continue our discussion on the dependence 
of Vt on distance, we first study the reproducibility of the 
measurements. We have repeated the experiments on dif- 
ferent contacts and for 3 samples. All data are plotted 
in Fig. [TJl. For all junctions Vt has a similar dependence 
on distance, and Vt has a maximum around 10"'^ Gq. 
However, there is a strong variation in the absolute value 
of Vt- For example around 10~^ Go, Vt scatters between 
1.4 and 2.2 V. Similarly, Vt usually differs for positive 
and negative voltage polarity in the same contact. This 
suggests that Vt is sensitive to the local shape of the elec- 
trodes. This is confirmed by closing the electrodes to a 
conductance > IOGq before starting a new measurement 
cycle. After this modification of the electrodes, Vt was 
found to have a different value. This is illustrated by the 
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FIG. 1: Electrical characteristics of clean gold junctions at 
Tf5i 5 K. a) Conductance versus width of the vacuum gap of 
a trained junction at a bias voltage oi V = 100 mV. From 
the exponential decay at large distance, an apparent barrier 
height is deduced of 4.2 ± 1 eV (Eq. [2]) Lia,liB]- b) Typical 
I-V curve in the tunneling regime. Here, Vt is -1.99 V for neg- 
ative voltage and 1.84 V for positive voltage (black dots), as 
obtained from: c) Fowler-Nordheim plot of the I-V character- 
istics for 34 different positions. After each curve the electrode 
separation is increased by 0.02 nm, resulting in a lower cur- 
rent (no offset is used). The black dots represent the minima, 
or K. Remarkably, Vt decreases with distance for wide tun- 
nel barriers while it increases for short tunnel barriers, d) Vt 
versus zero bias conductance for different contacts, measured 
on 3 different samples. Note the break in the scale between 
-1.2 V and +1.2V. The two curves marked by triangles and 
stars are measured on the same sample, but the latter was ob- 
tained after modifying the electrodes. The same holds for the 
two curves marked by squares and hexagons. For each con- 
tact, the apparent barrier height was measured and its value 
is given next to the data points (± 1 eV). 



filled triangles and filled stars in Fig. [TJi which represent 
two data sets that were obtained before and after such a 
reconfiguration. Actually, modifying the electrodes also 
changes the work function. This can be concluded by 
extracting the apparent barrier heights from the slopes 
of the exponential decay of the tunneling conductance 
for the individual contacts at large d (see Fig. [T^). The 
apparent barrier height is defined as 
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(f)^ = {hdlnG/dd)^/8m 



(2) 



For low voltages, (pa is expected to be almost equal 
to the work function of the electrodes [l^ [l^ . For our 
junctions (j)a was found to vary between 3.2 en 5.7 eV, 
which is in the range of typical values for gold electrodes 
with different crystal orientations [lO]. Remarkably, the 
contact corresponding to the green bullets in Fig. [TJi 
shows the highest values for both Vt and (pa- This brings 
us to the first conclusion: to obtain a reliable value of Vt 
the variation of the atomic configuration of the electrodes 
needs to be taken into account. 
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FIG. 2: Variation of the transition voltage with electrode 
distance. Filled symbols: measurements of Vt vs 1/d for 3 
different gold samples in vacuum. Dashed straight line: Vt 
vs 1/d as deduced from Eq. [l] (0 =5.4 eV). Dashed curved 
lines: Vt vs 1/d as calculated from the full Simmons model 
including the image forces. The strength of the image forces 
is varied by changing the prefactor ( in Eq. [1] from 0.4 to 
1. " x" and "+": data from large-area molecular junctions as 
reported by Beebe et al. "x": Alkane series, Cig-SH, Cig- 
SH, C12-SH, Cio-SH, Cs-SH, and Ce-SH. "+": phenyl series, 
TP-SH, BP-SH, and Ph-SH [3]. 

Let us now return to the distance dependence of T4 . In 
order to directly compare our measurements to the pre- 
dictions of Eq. [U we need to plot the data as a function of 
1/d. For this purpose, the origin in the position {d — 0) 
was defined by extrapolating the exponential part of Fig. 
[IJi (dashed line) to a conductance of 2e^//i. The crossing 
point is then set as the origin. As a result we obtain Fig. 
[21 There is a striking difference when comparing the ex- 
perimental data and the straight line expected from Eq. 
[TJ our data are not proportional to l/d. Instead, only 
a modest variation with d is found, with a maximum at 
« 3~^ nm~^. Clearly, the square barrier model (with a 
constant height 4>) does not give an accurate description 
of the data. 

As discussed by Simmons and applied by Huisman, a 
more realistic representation of the barrier takes image 



forces into account [tI, liclj . Image forces effectively round 
and lower the tunnel barrier. This can be described by 
the following equation for the shape of the barrier at a 
position X from the electrode surface: 



(x) = 00 - eV- - 0.58C- 
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where ^ is a variable we inserted to fit the strength of the 
image force [zlll^]- Actually, we note that Eq. [3] not only 
applies for image forces; in general, it may be regarded 
as a more elaborate description of the barrier. For ex- 
ample, one could also use a 1/r potential to describe the 
effective ion potentials of the outermost atoms in the two 
electrodes. At zero bias, the potentials of the opposing 
atoms will add up as —1 / x ~ 1 / {d — x)) to —d/{x{d~x)), 
resulting in the same type of expression as Eq. [3] 

Using the barriers described by Eq. [3l a different re- 
lation is obtained for Vt{d), which is plotted in Fig. [2] 
Interestingly, Eq. [3] yields a maximum in Vt at a cer- 
tain distance, in correspondence with the experimental 
data. We find the best qualitative agreement between 
the calculated values of Vt and the experimental values 
for a reduced value of the image force C = 0-4 — 0.7. 
This can be expected since the image forces for atomi- 
cally sharp electrodes are supposedly not as large as for 
infinitely large flat plates (as used by Simmons). Over- 
all, the rounded barrier gives a better description of the 
data. However, the variation of the experimental data 
with d is still much less pronounced than in the models. 

Let us first discuss the implications of our data for the 
interpretation of TVS. As described before, we want to 
make a comparison between the distance dependence of 
Vt for vacuum junctions and molecular junctions. For 
this purpose, we also include two data sets measured by 
Beebe et al. in Fig. [5] Js^. The upper data set (" x " ) , 
corresponding to alkanethiols of different lengths, has a 
negligible variation in Vt. This was ascribed to an almost 
constant HOMO-LUMO gap for different alkane lengths. 
The lower data set ("+") corresponds to (tt— conjugated) 
phenylene molecules. Compared to the alkanethiols, 
these molecules have a stronger dependence of Vt on 
d. This was attributed to a variation in the HOMO- 
LUMO gap, which is expected to decrease with increasing 
molecule length. Let us now compare the measurements. 
The two types of data (molecular junctions and vacuum 
junctions) have been measured at slightly different dis- 
tances. In contrast to our MCBJ measurements, the ex- 
periments of Beebe et al. were carried out on large-area 
junctions, in which many molecules are probed in paral- 
lel. In addition, the distance between the Fermi level and 
the nearest molecular level is lower than the work func- 
tion of the electrodes. As a result, the conductance is 
larger for the large-area molecular junctions which makes 
it possible to measure at larger distances or smaller val- 
ues of 1/d, respectively. Nevertheless, we find that the 
distance dependence of Vt for the molecular data does 
not differ significantly from that observed in our vacuum 
measurements. This is an important conclusion of this 
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paper. Taking into account (i) the measurement accuracy 
of the molecular data of approximately ±100 mV and (ii) 
the limited variation of Vt with d for the vacuum data it 
is not possible to distinguish molecular junctions from 
vacuum junctions just by measuring the distance depen- 
dence of Vt. However, considering the absolute values 
of Vt for conjugated molecules there is a clear difference 
with the vacuum data. For conjugated molecules, the re- 
ported values for Vt are much lower (0.6 V to 1 V) than 
the values found for the vacuum junctions (> 1.4 V) 0- 
This difference can thus be regarded as a signature of 
molecular junctions. 

Let us return to Fig. [2] and focus on the tunnel data 
only. Another way to interpret this graph, is that we have 
used TVS to test the tunnel barrier models described by 
Eqs. [I] and [31 In principle, any type of tunnel model 
could be tested, making TVS a potential tool for more 
detailed studies on electron transmission in a variety of 
tunnel systems. This is the third conclusion of our study. 
In our specific case, TVS allows us to conclude that the 
models in Eqs. [T] and |3] are incomplete, as they do not 
predict the observed weak Vt(d)-dependence. 

We note that other physical mechanisms affecting Vt 
are electrostatic forces [ly, [l^, [2^ and thermal ex- 
pansion. However, these can only lead to an additional 
lowering of Vt for very small gaps (< 0.25 nm). Hence, 
we should consider other mechanisms to explain the rel- 
atively weak variation of Vt with d. Actually, there is 
also another indication that the theory behind TVS is 
not yet complete. More complex TVS transitions have 
recently been reported in literature, where there is not a 
single well-defined minimum in the FN plot, but rather 
a plateau [l^, [IHl ■ These transitions have not been fully 
explained. In fact, we have also observed more complex 
transitions. In about 15% of the contacts, two minima 
appear close to each other (see supporting info). The 
fact that these complex structures appear for clean junc- 
tions without molecules indicates there are other physical 
mechanisms which could affect TVS. A possible explana- 
tion is due to the local density of states (LDOS) of the 



electrodes [2l[. If surface states are present, the LDOS of 
Au is strongly dependent on the electron energy. Thus, 
it is very well possible that a peak in the LDOS results 
in a minimum in the FN plot. Since the LDOS strongly 
depends on the crystal orientation and the exact shape 
of the electrodes, double minima may appear experimen- 
tally. Two minima generally only appear for one polarity 
of the voltage, indicating that the appearance of two min- 
ima is critically related to the exact properties of one of 
the two electrodes. 

Summarizing, we have performed TVS measurements 
on metal-vacuum-metal junctions. From the vacuum 
measurements we have learned the following aspects for 
TVS. First, Vt is found to be sensitive to the local shape 
of the electrodes. Thus, for a reliable value of Vt one 
should obtain information on the detailed configuration, 
or otherwise average over many configurations. Second, 
Vt is found to be sensitive to the distance between the 
electrodes. However, the relation between Vt and elec- 
trode separation d is much weaker than expected from 
standard models for vacuum tunneling and is compara- 
ble to the variation found for molecular junctions. Hence, 
it is not possible to use the variation of Vt on c? to dis- 
tinguish molecular junctions from vacuum junctions, as 
suggested recently, but the value of Vt provides a strong 
indication. Finally, we propose TVS as a simple but pow- 
erful test for models of vacuum tunneling. 
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